Numerical study of a three-state host-parasite system on the square lattice 
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We numerically study the phase diagram of a three-state host-parasite model on the square lattice 
motivated by population biology. The model is an extension of the contact process, and the three 
states correspond to an empty site, a host, and a parasite. We determine the phase diagram of the 
model by scaling analysis. In agreement with previous results, three phases are identified: the phase 
in which both hosts and parasites are extinct (So), the phase in which hosts survive but parasites are 
extinct (5*oi ), and the phase in which both hosts and parasites survive (Sow). We argue that both 
the 50-5*01 and S01-S012 boundaries belong to the directed percolation class. In this model, it has 
been suggested that an excessively large reproduction rate of parasites paradoxically extinguishes 
hosts and parasites and results in So- We show that this paradoxical extinction is a finite size effect; 
the corresponding parameter region is likely to disappear in the limit of infinite system size. 

PACS numbers: 87.23.Cc, 05.50.+q, 64.60.an 



I. INTRODUCTION 



In research fields ranging from ecology and epidemi- 
ology to sociology, it is important to clarify the effect 
of the interactions among species or phenotypes on the 
entire system. Stochastic interacting particle systems, 
in which each site on a graph takes either of the possible 
states and is flipped according to the states of other sites, 
are a useful tool for this purpose. A paradigmatic inter- 
acting particle system that describes disease spreading 
is the contact process (CP; also termed the susceptiblc- 
infected-susceptible model) [H-@|. 

Various interacting particle systems in complex net- 
works have been investigated recently [J, [|| . Neverthe- 
less, in an ecological context, organisms of different scales 
can be considered to live in a two-dimensional space, of- 
ten with a small interaction range. Therefore, it is in- 
structive to study models that are more complex than 
the CP on the Euclidean lattice 

0,0,0. A simple ex- 
tension of the CP in this direction is a three-state spatial 
host-parasite (HP) model that deals with an ecosystem 
comprising soil (empty sites), host species living on soil, 
and pathogen species (parasites) living on hosts. Phase 
transitions and oscillations in similar models have been 
studied from the perspective of statistical physics [81412} . 



'Electronic address: hascgawa® st at . t . u- toky o. ac . j p | 

t Electronic address: konno@y nu.ac.jp| 

t Electronic address: masuda@mist.i.u-tokyo. ac.jp 



Sato et al. [13( analyzed the HP model on a square 
lattice. They showed by means of the improved-pair ap- 
proximation (i-PA) and numerical simulations that a very 
high reproduction rate of parasites results in the extinc- 
tion of both hosts and parasites. This phenomenon is 
called parasite-driven extinction [l3l [l4|. An intuitive 
explanation for this paradoxical behavior is that para- 
sites replace hosts so quickly that hosts get extinct, which 
eventually results in the extinction of parasites. A similar 
paradoxical behavior, i.e., a decrease in the number of a 
species caused by an increase in its fertility, is observed in 
other models, where a sort of rock-scissors-paper compe- 
tition is prevalent among three species [15l423| . However, 
the current understanding of the phase diagram of the 
HP model is not comprehensive, because parasite-driven 
extinction cannot be predicted by mean field approxima- 
tion and pair approximation (PA) 0, HH . 

In this paper, we numerically investigate the phase di- 
agram of the HP model on the square lattice. In par- 
ticular, we use large lattices and investigate the effect of 
the system size on parasite-driven extinction. The ob- 
tained phase diagram is shown in Fig. Q] We argue that 
two transition boundaries (solid lines in Fig. [T]) belong to 
the directed percolation (DP) universality class. Another 
transition boundary (dotted lines) is not characterized by 
the DP universality class, and its location depends on the 
system size. We claim that the parasite-driven extinction 
phase is a finite size effect and that the phase diagram is 
qualitatively the same as that obtained by the PA rather 
than that obtained by the i-PA. 
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FIG. 2: Transition rules of the HP model. Solid and dashed 
lines represent deaths and births, respectively. The values 
indicate the transmission rates, and m denotes the number of 
neighbors of a site in state i. 



FIG. 1: Phase diagram of the HP model. The solid blue line 
represents the boundary between So and Soi (i.e., Aoi = A§i). 
The solid red line represents the boundary between Soi and 
£012 an d is drawn on the basis of the data shown in Table [T] 
Above the dashed lines, which correspond to L — 300, 500, 
and 1000 from the bottom to the top, the parasite-driven 
extinction occurs frequently. 

II. MODEL 

The HP model on the square lattice Z 2 is defined 
as a continuous-time Markov process with state space 
{0, l,2} z2 (H H HiHiJ. Each site ta kcs one of the 
three states 0,1, and 2, which represent an empty site, a 
host, and a parasite, respectively. The rules for the state 
transition are depicted in Fig. [5] A host and a parasite 
die at rates d\ and cfe, respectively. For simplicity, we 
set d\ = c?2 = 1- The occurrence of death at any site 
is independent of the states of the neighboring sites. In 
contrast, reproduction of hosts and parasites depends on 
the states of the neighbors. A host emerges at an empty 
site i at rate Aoini(i), where nx(i) is the number of hosts 
in the neighborhood of site i. A host at site i turns into 
a parasite at rate Ai2?T.2(«), where 112(1) denotes the num- 
ber of parasites in the neighborhood of site i. We vary 
the values of Aoi and A12 in the numerical simulations. 
Because parasites feed on hosts, the HP model allows the 
following three phases in the stationary state: 

(i) phase So, in which hosts and parasites are extinct, 

(ii) phase 5*01, in which hosts survive and parasites are 

extinct, and 

(iii) phase So 12, m which both hosts and parasites sur- 
vive. 

The HP model with A12 = is equivalent to the CP. 
In the CP, each site takes either state or 1, and a death 
event (1 — >• 0) and a reproduction event (0 — > 1) at site 
i occur at rate di = 1 and Aoini(i), respectively. In 
the case of the CP on the square lattice, So and Soi are 
realized when Aoi is respectively smaller and larger than 
Agj « 0.4122 [l|. 



The phase diagram of the HP model on the square 
lattice has been examined using the mean field approx- 
imation fl3|; th e PA, which accounts for pairwise state 
correlation [2J|; and the i-PA, which calibrates the PA 
to account for the aggregation of the same species in the 
space [l3|, [lij . All of the three approximations predict 
the existence of the three phases of the model. In the 
mean field approximation and the PA, the system is in 
So if Aoi is less than a critical value that is independent 
of Ai2- Otherwise, the system is in Soi (S012) when the 
value of Aoi and A12 is sufficiently small (large). In the 
mean field approximation and the PA, the boundary be- 
tween Soi and S012 is given by A12 = Aoi/(4Aoi — 1) 
and A12 = (12(A i) 2 + 4A i)/(36(A i) 2 - 4A i - 3), re- 
spectively [2J]. In particular, only So and S012 exists 
when A12 — > 00 in the mean field approximation. In the 
PA, when A12 — > 00, So, Soi, and S012 appear in this 
order in the PA as Aoi increases. The phase diagram ob- 
tained from the i-PA is qualitatively distinct from those 
obtained from the mean field approximation and the PA. 
When A12 is large, the i-PA predicts So regardless of the 
value of Aoi- This result corresponds to the numerical 
observation that a large reproduction rate of parasites 
induces extinction of hosts and parasites [HI, flij ]. We 
call this phenomenon the parasite-driven extinction. The 
mean field approximation and the PA do not predict the 
existence of the parasite-driven extinction. 



III. DP TRANSITION ON THE S01-S012 
BOUNDARY FOR SMALL X 12 

In this section, we numerically examine the boundary 
between Soi and S012 for small values of A12 (the red solid 
line in Fig.[T]). We carry out Monte Carlo simulations for 
the HP model on the square lattice with N = L x L 
sites, where L = 300. Periodic boundary conditions are 
assumed. We run 500 realizations for fixed Aoi and A12. 
At the beginning of each realization, each site indepen- 
dently takes state 0, 1, or 2 with equal probability. We 
adopt an event-driven update algorithm in which we se- 
lect one out of all the possible events to occur with the 
appropriate probability for each time step. Then, we in- 
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crement the time by an appropriate amount. 

First, we focus on the limit Aoi — > oo, where an empty 
site adjacent to a host is instantaneously replaced by the 
host. A cluster of empty sites survives only when they 
are surrounded by a shell of parasites. When A12 is small, 
parasites rarely form such a shell. Then, the HP model 
behaves like the CP, where empty sites and hosts in the 
HP model collectively correspond to the susceptible sites 
(i.e., state 0) in the CP. Because many spatial stochas- 
tic processes including the CP undergo a DP-typc phase 
transition d H^-ii], we expect that the HP model also 
undergoes a DP-type transition from S01 to 1S012 as A12 
is increased to cross « Af^ ~ 0.4122. The time courses 
of the mean density of parasites (p2)(~t) are shown in 
Fig. [3ja) for various values of A12, where (•) denotes the 
average over all the realizations. At A12 = Aj 2 w 0.4129, 
we obtain 

(p 2 )(t)Kt- s . (1) 

From Fig. EJb), which shows the plotting of the local 
slopes of {p2){t) 1 we obtain 8 « \og{p2)(t)/\ogt « 0.451, 
a value indicative of the DP univers ality class We 
also derive 6 via dynamic scaling [27i [Hf, i.e., by fitting 
the following scaling form: 



(P2){t)~t 



where 



AA 



12 



Pi AX 12 t 



A12 — A, 



(2) 



(3) 



The critical exponent <5 is given by 5 — fijv\\. The re- 
sults of the dynamic scaling with the known critical ex- 
ponents for the (2+l)-dimensional DP universality class 
/? « 0.583 and v\\ w 1.295 Q are shown in Fig.^c). The 
data for different values of A12 collapse onto a single curve 
separately for subthreshold and suprathrcshold values of 
A12. This result also supports that the transition belongs 
to the DP universality class. 

If Aoi is finite and sufficiently large, we can numerically 
obtain the transition points and the critical exponents in 
the same manner. On the critical line, {p2)(t) shows a 
power law decay with t, as shown in Fig. Eta). When 
Aoi ^ 0.68, the dynamic scaling yields the DP critical 
exponents at each examined transition point. The lo- 
cations of several points on the <S'oi-5'oi2 boundary are 
shown in Fig. [T] and Table Q] 

We postpone the analysis of the case Aoi < 0.68 to 
Sec. ED 



IV. DEPENDENCE OF BOUNDARY BETWEEN 
5012 AND THE PARASITE-DRIVEN 
EXTINCTION REGION ON A 12 

Parasite-driven extinction may occur for large A12 
[l3l [l4| . Figure Hl^b) shows the surviving probability of 
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FIG. 3: DP phase transition at Aoi — > 00 and A12 ~ A^2- 
(a) Time courses of (p2){t) and (b) local slope S of (p2)(t). 
The different lines from the top to the bottom correspond 
to A12 = 0.4079, 0.4089, . . ., and 0.4179. (c) Dynamic scaling 
(Eq. ((2}) for the data shown in (b). 



hosts Pi (t) and that of parasites P2 (t) for some large val- 
ues of A12 and fixed values of Aoi = 10 and L = 300. If 
Pi (t) approaches zero rapidly, the parasite-driven extinc- 
tion is considered to have occurred. If the transition from 
So 12 to the parasite-driven extinction belongs to the DP 
universality class, Pi(t) or Pzit) should decay geometri- 
cally on the phase boundary and exponentially for A12 
slightly larger than the critical value. 

However, Fig. H^b) indicates that this is not the case. 
Whether extinction of hosts and parasites occurs or not 
is determined at an early stage, where hosts are rapidly 
replaced by parasites, resulting in a rapid decrease in the 
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TABLE I: Several points on the S01-S012 boundary. 



A01 


0.509 0.543 


0.591 
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FIG. 4: (a) Time courses of (p 2 )(£) with A01 = 10. The 
lines from the top to the bottom correspond to A12 = 
0.430, 0.432, . . ., and 0.450. (b) Surviving probability of hosts 
Pi (t) (dashed lines) and that of parasites P2 (i) (solid lines) 
with A01 = 10. The lines from the top to the bottom corre- 
spond to A12 = 5.9, 6.3, 6.7, 7.1, and 7.5. We set L = 300 in 
both (a) and (b). 



number of hosts. If the hosts die out, the parasite-driven 
extinction takes place. In contrast, if hosts survive the 
initial stage, which occurs with a low probability, the 
hosts recover from near extinction. In this case, hosts 
and parasites are likely to coexist for long time. The 
value of A12 affects the probability that the hosts survive 
rather than the rates at which the number of hosts and 
parasites decay. 

We state that the parasite-driven extinction is a finite 
size effect. In order to confirm this statement, we mea- 
sure the probability of the parasite-driven extinction as a 
function of linear lattice size L. Because the transient is 
short, as shown in Fig. E^b), we measure the fraction of 
realizations among 2000 realizations in which both hosts 
and parasites are extinct at t = 100. Figure EJa) shows 
the extinction probability for a range of values of A12 at 
A01 = 10 and L = 100, 200, 300, 500, 700, and 1000. The 
extinction probability indefinitely decreases with L. The 
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FIG. 5: (a) Relationship between the extinction probability 
and A12 when L = 100, 200, 300, 500, 700, and 1000 (from left 
to right), (b) Dependence of A{ 2 (Aoi, L) on L. We set A01 = 
10 in both (a) and (b). The number of realizations for a given 
combination of A12 and L is equal to 2000. 

value of A12 that has an extinction probability of 1/2, 
denoted by A12 = A{ 2 (Aoi,£), is plotted against L in 
Fig. E{b). It is observed that A{ 2 (Aoi,£) oc InL. Loga- 
rithmic scaling is also observed at other values of A01. In 
Fig. [1] we show A{ 2 (Aoi, L) for some values of A01 and L 
(dotted lines). 

The results obtained in this section indicate that the 
parameter region of parasite-driven extinction indefi- 
nitely shrinks as L increases. This system-size depen- 
dence is distinct from the dependence of the critical value 
on L in the usual phase transitions, which is convergent 
in the limit of infinite system size. 



V. DP TRANSITION ON THE S -S i 
BOUNDARY IN THE LIMIT A12 -> 00 

When A12 is sufficiently large, the mean field approxi- 
mation predicts that the system transits from So to S012 
as A01 increases [l3[. The PA predicts that the system 
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transits from So to Soi and then to 6*012 as A01 increases 
[23 | . The i-PA predicts that the system is in Sq irrespec- 
tive of the value of A01 (see Figure 1 in [Hj])- To analyze 
this apparent contradiction, we carry out simulations in 
the limit A12 — > 00. 

Irrespective of the value of A12, it seems that A01 must 
be larger than A^ for hosts to survive. Therefore, we 
start by examining the case A01 ~ A^. When A01 ~ Aqi 
and A12 — > 00, a host adjacent to a parasite is instan- 
taneously invaded by the parasite. In such a case, if 
we start numerical simulations on the equal fraction of 
empty sites, hosts, and parasites, then the number of 
hosts, if they survive at all, becomes small at the very 
beginning of a run. For example, the averaged number 
of hosts on the 300 x 300 square lattice decreases from 
30000 to w 60 after a short time. It may not be suitable 
to measure the decay of the expected number of hosts, 
which would be (pi)(t) oc t~ s on the critical line; this is 
because such a measurement necessitates the existence of 
a sufficiently large number of hosts at the beginning of a 
run. 

Another numerical method for estimating the transi- 
tion point and critical exponents is to measure the time 
courses of the system starting from an almost absorbing 
configuration [301 ] . For example, we observe the power 
law behavior of the surviving probability, the number of 
active sites, and the mean spreading at the transition 
point, if we run the CP starting from a single active site. 
Therefore, we assume that the initial configuration of the 
HP model contains just one host. The other sites are ei- 
ther empty or parasites with a probability of 0.5. 

With this one-host configuration, the mean number of 
hosts follows the power law {N\){t) oc t e at A01 ~ A^, as 
shown by the solid lines in Fig. Ula). On the other hand, 
parasites rapidly become extinct (dashed line). The sur- 
viving probability of hosts also follows the approximate 
power law P 1 (t) oc t~ s ' in the same parameter range 

(Fig. EDO) 

At A01 r 
represented by 



Aq 1; we adopt the dynamic scaling ansatz [l[ 
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where 



AAqi — A01 — Aqi- 



(6) 



This dynamic scaling ansatz explains the data shown in 
Figs. EI a ) and Eljb), respectively. The fitting results 
with the DP exponents 6 re 0.229 and 5' = 8 0.451 Q 
(Figs. Etc) and Eld)) suggest that the transition from So 
to Sqi at A01 = \qi and A12 — ► 00 is of the DP type. We 
consider that this phase transition is independent of the 
value of A12. This result qualitatively agrees with that 



FIG. 6: Transitions at A01 ~ A01 and A12 — > 00. (a) Time 
courses of (Ni)(t) (solid lines) and (N2)(t) (dashed lines that 
almost overlap each other), (b) Surviving probability of hosts 
Pi(t). (c) Dynamic scaling (Eq. @) for the data shown in 
(a), (d) Dynamic scaling (Eq. (JSJl) for the data shown in (b). 
The lines correspond to A i = 0.4082, 0.4092, . . ., and 0.4182 
from the bottom to the top. The number of realizations for a 
given A01 is equal to 10 7 . 
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FIG. 7: (a) Time courses of (Ni}(t) when A i « A§j. (b) 
Surviving probability of hosts Pi(t). (c) Dynamic scaling 
(Eq. ((4))) for the data shown in (a), (d) Dynamic scaling 
with Eqs. (J5]) and ((Sj applied to the data shown in (b). The 
lines correspond to A i = 0.40821,0.40921,..., and 0.41821 
from the bottom to the top. The number of realizations for a 
given Aoi is equal to 10 7 . 



obtained from the PA but not that obtained from the 
i-PA. 

With the random initial configuration, we observe 
(Ni)(t) and Pi{t) instead of (pi)(t) and obtain the same 
results as those shown in Fig. [6] (Ni)(t) and Pi(t) decay 
geometrically at Aoi ~ Ag l5 as shown in Fig. [T^a) and 
E^b), respectively. The dynamic scaling (Eq. with 
the DP exponents fits (Ni)(t) shown in Fig. [TJa) well 
(Fig. [TJc)). On the other hand, dynamic scaling of P\(t) 
(Eq. (O) fails because the number of surviving hosts af- 
ter a short time is greater than one. To circumvent this 
case, we assume that the surviving hosts are located away 
from each other and grow independently on the lattice. 
We denote the surviving probability of a specified host 
by P^ mglc (i). Then, we approximate P\{t) as 

a(t)«l-(l-if n * U W l , (7) 

that is, 

Pf nglo (t) ^-(l-Pi^)) 1 /", (8) 

where n is the mean number of surviving hosts after a 
short time. By replacing Pi(t) in Eq. (|5|) by Pf lngle (i) and 
using the DP critical exponents, we obtain a reasonable 
scaling, as shown in Fig. EJd). 

VI. S012 PHASE IN THE LIMIT A 12 —¥ oo 

When A12 — > oo, either the random initial configura- 
tion or the one- host configuration yields So or Soi, but 
not 5*0121 for any value of Aoi- This remains the case for 
at least up to L = 1000. The apparent absence of S012 
may be because there are initially too many parasites. 
In the case of a large A12, parasites replace hosts in a 
short time, which is likely to lead to the extinction of the 
parasite. 

To examine the possibility of S'012 at A12 — > 00, we 
adopt the one-parasite configuration, where the remain- 
ing sites are either empty or occupied by the host with 
a probability of 0.5. With this initial configuration, we 
find that both hosts and parasites can survive when L is 
large and Aoi is within a certain range. When L < 400, 
neither hosts nor parasites survive. 

Time courses of the number of parasites are shown in 
Fig. ||for L = 700 and three values of Aoi- As Aoi in- 
creases within this range, the basal number of parasites 
in a short run increases, but the amplitude of the damped 
oscillation in the number of parasites also increases. If 
Aoi is sufficiently large, the amplitude of the oscillation 
is so large that the parasites are likely to disappear in 
the first cycle of the oscillation (Fig. HJc)), whereas the 
basal number of parasites is larger than that in the case 
of a smaller Aoi (e.g., Fig.[5ta)). We remark that, for re- 
lated spatial stochastic processes, sustainable oscillations 
[3TI . |32| and absorption to the unanimity state owing to 
the blowing out of oscillations [33| were reported as finite 
size effects. 
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FIG. 8: Time courses of N 2 (t) when A i = (a) 0.515, (b) 
0.530, and (c) 0.545. We set A12 oo and L = 700. Each 
colored line represents a single run, and the results for 4 runs 
are shown in each panel. 



The stationary density of the parasites averaged over 
the surviving runs, denoted by (/92}surv, is shown for some 
large values of L in Fig. HJa). Here (.} SU rv indicates the 
average over realizations in which parasites survive after 
a transient of length 1500. We observe that (p2)surv is 
positive for Aoi > 0.509 and converges to a certain value 
for Aoi > 0.518. We did not determine the transition 
point and the critical exponents by a scaling argument 
for (/?2)surv m terms of Aoi because {p2)smv is too small 
for Aoi ~ 0.509. To support the existence of the So 12 
phase in the limit L — > 00, we measure the fraction of 
surviving runs for various system sizes. As shown in 
Fig. Elb), the fraction of surviving runs increases with 
L for Aoi > 0.509. This result supports the fact that 



(a) 

<P2> S1 

0.008 
0.006 
0.004 
0.002 
0.000 

(b) 

% surviving runs 
l.Or 

0.8 

0.6 

0.4 

0.2 

0.0 



0.510 0.515 0.520 0.525 0.530 0.535 



0.51 0.52 0.53 0.54 



0.55 



FIG. 9: (a) Stationary parasite density (p2)surv averaged over 
the surviving runs in the limit A12 — ¥ 00. (b) Fraction of the 
surviving runs. We set L = 500 (triangles), 600 (diamonds; 
only in (b)), 700 (squares), and 900 (circles). The number of 
realizations for a given combination of Aoi and L is equal to 
100. 



S012 exists for Aoi > 0.509 in the limit L — > 00. As Aoi 
increases even further (i.e., Aoi > 0.524), the fraction of 
surviving runs decreases. The parasite-driven extinction 
for a finite system size gets eminent in this range of Aoi; 
this parasite-driven extinction is caused by the increasing 
magnitude of damped oscillations. Similar to the results 
shown in Sec. IIV1 the parameter region for the parasite- 
driven extinction depends on the system size and is likely 
to disappear in the limit L^oo. We also observed that 
the results in the case of finite A12 > 2 are qualitatively 
the same as those in the case of A12 = 00. 

Finally, we examine the S01-S012 transition line for 
large A12. In this case, we do not obtain a data collapse 
by the dynamic scaling based on the relaxation of the 
system, as shown in Fig. [TUt a) for A12 = 4. Therefore, 
we attempt the dynamic scaling for the parasites in the 
manner similar to that employed in Sec.fV] Consider the 
neighborhood of the S01-S012 transition point for a large 
fixed Ai2- With the one-parasite configuration, a parasite 
would quickly invade hosts at an early stage. In this case, 
the growth rate of the parasite is fairly insensitive to Aoi- 
Therefore, the scaling argument would not apply. 

To avoid such an initial growth of parasites and obtain 
a clear scaling of (A^Ht), we proceed as follows. First, 
we start a simulation from a mixture of independently 
distributed empty sites and hosts with the equal density 



8 



(a) 

<iv,>(t), <JV 2 >(t) 




10 10 2 10 : 



(b) 

<iV 2 )(t) 




FIG. 10: (a) Time courses of (Ni)(t) (dashed lines) and 
(N2)(t) (solid lines) for the random initial configuration. 
The lines from the bottom to the top correspond to Aoi = 
0.58,0.582,..., and 0.60. The number of realizations for a 
given Aoi is equal to 50000. (b) Time courses of (N2)(t) for 
the modified initial configuration. The lines correspond to 
Aoi = 0.57, 0.575, . . ., and 0.615 from the bottom to the top. 
The number of realizations for a given Aoi is equal to 20000. 
We set A12 = 4 and L — 300 in both (a) and (b). (c) Dynamic 
scaling for the data shown in (b). As the scaling function, we 
use Eq. © with (Ni)(t) replaced by {N 2 )(t). 



(i.e., 0.5 each). After the system has approached a steady 
Soi state, we replace a randomly chosen empty site with 
a parasite and continue the simulation until the station- 
ary state is reached. Figure HUT b) shows the time course 
of (N2)(t) for A12 = 4 and various values of Aoi, where 
the single parasite is added at t = 0. Near the transition 
point, Aoi ~ 0.591, (A?2)(i) seems to follow a power-law. 
The data for different values of A12 collapse onto a single 
curve with the DP critical exponents, separately for sub- 
threshold and suprathreshold values of A i (Fig. ITUtc)). 
Figure [TUtc) suggests that the transition belongs to the 
DP universality class. 

Note that (JVa)(i) above the transition point saturates 
owing to a finite size effect. It is difficult to deter- 
mine critical properties for large values of A12 because 
we would need increase L to perform the dynamic scal- 
ing. Nevertheless, we believe that the Soi-Soi^ transition 
belongs to the DP universality class even for larger A12. 



VII. SUMMARY 

We carried out numerical simulations for a three-state 
host-parasite model on the square lattice. The obtained 
phase diagram is shown in Fig. [T] Our numerical re- 
sults suggest that the Sq—Sqi boundary and the Soi - £012 
boundary are of the DP universality class. The parasite- 
driven extinction occurs for large Aoi and large A12 in 
relatively small systems. However, for a sufficiently large 
system, the three states coexist in the parameter region 
where the parasite-driven extinction occurs for a small 
system. Therefore, the parasite-driven extinction is a fi- 
nite size effect. This prediction is consistent with the 
phase diagram obtained from the PA but not with phase 
diagrams obtained from the mean field approximation 
and the i-PA. 
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